Time-dependent density-matrix renormalization-group: A systematic method for the 
study of quantum many-body systems out-of-equilibrium 
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The density-matrix renormalization-group (DMRG) algorithm is extended to treat time- 
dependent problems. The method provides a systematic and robust tool to explore out-of- 
equilibrium phenomena in quantum many-body systems. We illustrate the method by showing 
that attractive interactions enhance the tunneling current between two Luttinger liquids, whereas 
repulsive interactions suppress it, in qualitative agreement with analytical predictions. Enhancement 
of the transport current through a quantum dot in the Kondo regime is also exhibited. 
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The study of out-of-equilibrium phenomena in quan- 
tum many-body systems has mainly focused on steady 
state properties, such as the conductance of quantum 
dot systems p], p| and of point contacts between two Lut- 
tinger liquids g7 Hj- More recently, attention has been 
paid to time-dependent phenomena in these systems]^, 



|, g, |, |, [TO], [TTJ, [L2| |TJ, |TJ. One goal of the work 
has been to identify the hierarchy of time scales which 
arise when quantum-many body systems are driven out 
of equilibrium by a variety of time-dependent perturba- 
tions. This information is relevant for experiments in- 
volving potentially important technological applications 
of quantum dots and junctions built from carbon nan- 
otubes, and quantum wires which behave as Luttinger 
liquids (LL)g, However, many studies have em- 

ployed either uncontrolled approximations, or exact solu- 
tions at special values of the model parameters which are 
rather remote from experimentally accessible systems. 

By contrast, the density-matrix renormalization-group 
(DMRG) algorithm [DJ has been successfully em- 
ployed to calculate many static and frequency-dependent 
properties of quasi one-dimensional quantum many-body 



systems) 17 . The numerical method does not depend on 
whether the system is integrable or not, and in many 
cases it provides essentially exact results as errors in- 
duced by truncating the Hilbert space are controlled in 
a systematic way. 

We employ the so-called infinite-size dmrg (T6[ [n) as 
it treats the center of the chain, the region where we fo- 
cus our attention, most accurately. Our starting point is 
a quantum Hamiltonian representing a chain of L sites. 
If we denote by D the dimension of the Hilbert space on 
each site, then the dimension of the Hilbert space Tt is 
D L . The DMRG algorithm organizes the Hilbert space 
into blocks, Bl • »Br, where represents the Hilbert 
space of the two central sites (of dimension D) and Bl, 
Br are respectively the blocks representing the Hilbert 



space of the remaining left and right sites, each of di- 
mension £) L / 2 ~ 1 . An initial small chain of L = 4 sites is 
first enlarged to 6 sites by cutting it in half and inserting 
two sites at the center. As this process is repeated, the 
Hilbert space grows exponentially and eventually exceeds 
a limit, (D * M) 2 , beyond which it is truncated. As the 
calculation may be systematically improved by increas- 
ing the size of the retained block Hilbert space, M, up 
to limits imposed by computer memory and speed, it is 
straightforward to check for convergence in the observ- 
ables. The end result of the algorithm is an approxima- 
tion to the exact many-body ground state I^Jf 11110 ), the 
ground state energy Eq, and a truncated representation 
of the Hamiltonian, i?trunc- 

The extension to time-dependent problems, which 
we call time-dependent DMRG (TdDMRG), takes the 
infinite-size DMRG algorithm described above as its 
starting point. After the penultimate iteration, when the 
chain has L — 2 total sites, either a quantum dot paired 
with an ordinary site, or two sites linked by a tunnel- 
ing junction (depending on the choice of problem to be 
studied) are inserted at the center of the chain, yielding 
a chain with a total of L sites. Finally, a time-dependent 
perturbation, H'(t) is added to H tTU nc, and the resulting 
time-dependent Schrodinger equation is then integrated 
forward in time: 

= [(tftruac - E ) + H'{t)\ |*(t)) . (1) 

The initial state is chosen to be the ground state of 
the unpertured truncated Hamiltonian, \^(t = 0)) = 

To illustrate the method, we apply the TdDMRG al- 
gorithm to two kinds of systems: quantum dots and tun- 
neling junctions between two chains of one-dimensional 
interacting fermions. We first consider the easier problem 
of spinless fermions (with small on-site Hilbert space di- 
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FIG. 1: Transport current J(t) for the spinless quantum dot 
system denned by Eq.(^). The plots (which are offset in 
the vertical direction for clarity) illustrate systematic conver- 
gence towards the exact independent-particle solution with 
increasing block dimension M. Time is measured in units of 
r m = h/w and J(t) is expressed in units of e/r w . The param- 
eters are t q = 0.15w, e q = — 0.25«j, and L = 64 sites. At time 
t — 4t w the chemical potential of the leads is shifted by a 
bias of Sfi — —0.25w within a rise time of t 3 = 0.1 t w . The 
period of the ringing oscillations is given by the expression 
reported in Ref. 0: Ata = 2ttH/\S^r — e q \ = 4ir r w . 



mension D — 2) before turning to a more realistic model 
involving electrons with spin. 

1. Spinless Fermions: We turn first to the prob- 
lem of a quantum dot coupled to two leads. A similar 
problem has been studied by Wingreen et al.(|] in the 
thermodynamic L — > oo limit. In terms of the fermion 



and the num- 



creation and annihilation operators cj , , 

ber operators rij = c^Cj, the time-independent part of 
the Hamiltonian for the quantum dot system reads: 



H, 



qdot 



w 
~2 



E 



[ c ]+i c j 



+ H.c. 



H.c. 



(2) 



In this equation q denotes the location of the right central 
site where the quantum dot is located, q = L/2 + 1. Here 
t q is the hopping amplitude between the dot and either 
lead, and w > is the half band-width of the leads, 
which arc half-filled. The local chemical potential at the 
quantum dot site, e q , is generally nonzero. It may be 
viewed as the energy of a localized level. 

The time-dependent perturbation is a differential bias 
applied between the left and right leads: 



H'(t) = -6fi R (t)N R - d/i L (t)N L 



(3) 



We typically set J^i(t) = —Sfin(t) = Sfi 8(t — t ), where 
9(t) is a smoothed step function with rise time t s : 6{t) = 
l/[exp(-t/t s ) + 1]. Operators N R = Xf 



and A^l 

(which has a similar expression) count the number of 
fermions in the right and left leads. 

It is convenient to gauge-transform the bias into a lo- 
cal, but time-dependent, hopping amplitude by means 
of the transformation |t(t)) = e i*(t)[JVR-JV L ]/^ |*(t)), 

where &(t) — Sfj,L{t') dt' . Thus the real-valued hop- 
ping amplitude t q is replaced by a complex-valued one: 
t q (t) — t q e^^^ h , and the differential bias, which acts on 
all the sites in the blocks, is absorbed into an operator 
that acts only locally on the central region of the chain 
where the infinite-size DMRG is most accurate. 

The expectation values of the currents along the cen- 
tral links may be evaluated during the course of the wave- 
function's time-evolution by calculating: 

J q , q -i(t) = - fuel* t q (t) {*{t)\c\c q ^{t))} . (4) 

Following Refs. || and || we introduce an average trans- 
port current defined as: J(t) = \[J q . q -i{t) + Jq+i, q {i)]. 
The results of our calculation are shown in Fig. |l| The 
exact solution for non-interacting systems is obtained 
by integrating the equations of motion for the quanti- 
ties jij = (c\(t)cj(i)) forward in time for systems of the 
same finite size L; therefore the results can be compared 
directly with the TdDMRG calculations. 

Another problem of recent interest is transport 
through a point contact or junction between two Lut- 
tinger liquids. As the DMRG algorithm works equally 
well for interacting leads, we may turn on nearest- 
neighbor interactions VrijUj+i between the spinless 
fermions (again q = L/2 + 1): 



H 



junct 



-f E [4+1^- +h. 



tq [cjc 9 _i + H.c] 



(5) 



At half-filling the leads are metallic for |V| < w and ex- 
hibit LL behavior ]l9|; otherwise a charge-density wave 
(CDW) forms and an insulating gap opens up. Kane and 
Fisher || studied this model in the continuum limit for 
V < w. Within perturbative RG and bosonization they 
concluded that intralead interactions strongly affect the 
transmittance of the junction. For attractive interactions 
(V < 0), the transmittance is enhanced whereas for re- 
pulsive potentials (V > 0) it is suppressed. 

Results for the two systems of spinless fermions are 
shown in Figs, [j] and ^. In both cases the short time 
behavior of the current compares well with the available 
exact independent-particle results, even for rather small 
block Hilbert space dimension M (see also Fig. || below) . 
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FIG. 2: Current across a junction between two systems of one- 
dimensional interacting electrons in response to a step bias 
applied between the two leads. Interactions strongly affect 
the currents, in qualitative agreement with the predictions of 
Ref. ^[ The junction consists of a weakened link with t q — 
0.125 w lying at the center of a chain of L = 64 sites. At 
time to — 5 t w the bias is shifted by 8fi = 0.0625 w in a 
rise time of t s = 0.1 t w . The label "linear response" refers to 
the application of a much smaller bias Sfi = 6.25 x 10 -4 w. 
Lines are for the more accurate M = 512 truncation (or, in 
the non-interacting V — case, depict the exact solution) 
while data points are for M = 256. For clarity, currents have 
been multiplied by factors of 2, 20, and 800 as indicated. 
Note the good convergence with increasing blocksize M, and 
with the exact result in the non-interacting case. The units 
of current and time are the same as in Fig. [j]. The reflection 
time t la fi = 64t w is longer than the time period shown. 



This may seem surprising, as only the ground state was 
targeted in the DMRG portion of the algorithm. How- 
ever it is important to note that the DMRG algorithm is 
not a true renormalization-group method as the reduced 
density-matrix, rather than energetics, determines which 
states are projected-out of the Hilbert space. Over short 
times the perturbation can only create a limited number 
of excitations about the ground state, so the time-evolved 
state remains closely connected with the ground state 
and can be accurately represented with a superblock of 
relatively small Hilbert space dimension. Fig. [I] also illus- 
trates the systematic convergence with increasing block- 
size M. By retaining more Hilbert space we are able 
to extend the accurate solution forward in time. Never- 
theless even short-time behavior often suffices to reveal 
interesting physics. For example the initial behavior of 
the junction current depicted in Fig. ^| is consistent with 
the prediction of Kane and Fisher ||: An attractive in- 
teraction enhances the current J(t) whereas a repulsive 
interaction suppresses it. It is interesting to observe that 
in the insulating CDW regime V > w current oscillations 
develop after the initial transient. We emphasize that the 



TdDMRG method enables us to access short-time out-of- 
equilbrium phenomena in an accurate way, and thus goes 
well beyond bosonization. 

Several other aspects of the TdDMRG method are 
worth mentioning. First, as the numerical method nec- 
essarily addresses only systems of finite size, the current 
averaged over long times must always vanish, regardless 
of the bias. This is so because any current pulse created 
by the bias eventually reaches the open ends of the chain 
within a time of order i re fl = La/vp, where vp = w a/h 
is the Fermi velocity and a the lattice spacing. At that 
point the pulse is reflected back, leading to oscillations 
in the current. To be useful, model parameters must be 
chosen such that the physics of interest occurs over time 
scales shorter than i rc fl. There is no limit on the size 
of the applied bias: large biases can be applied, driv- 
ing the system far outside of the linear-response regime, 
though we find that larger biases generally require larger 
blocksizes M to maintain the accuracy of the solution. 
Sizeable biases were applied in the above two model prob- 
lems, Figs. |]and The non-linear effects of the bias are 
illustrated by the response in the V > w regime shown 
in Fig. ||, as the form of the oscillations depends strongly 
on the size of the applied bias. 

2. Spinning Electrons: To illustrate the applica- 
tion of the TdDMRG method to a more realistic problem 
involving spinning electrons, we turn now to the study 
of the one-impurity Anderson model described by the 
Hamiltonian: 



H, 



qdot 



W 

~2 



E 

+ e q n q + U n^ q n lq 



t<7 
C j+l C j<? 



H.c. 



(6) 



Here a =t,| labels the spin and n„ 



the total number of electrons at the dot site. The one- 
impurity Anderson model is a minimal model to de- 
scribe a quantum dot with spinning electrons coupled 
to leads B. The model also arises in studies of strongly 
correlated electrons systems in the limit of infinite spa- 
tial dimension, in particular within dynamical mean field 
theory (DMFT)Q. The method presented here can be 
adapted to the calculation of the imaginary-time Green's 
functions needed in the DMFT calculations |H| . 

In this case only the bias of the left lead is shifted, 
H'{t) = —Sfj,L{t)Ni,, but again this shift can be gauge- 
transformed into a time-dependent hopping amplitude. 
Results for the initial transport current J(t) are shown 
in Fig. |. Here e q = -0.25w and the width A = 0.125w; 
thus in the large-U limit the quantum dot is in the Kondo 
regime. For these parameters the Kondo scale is esti- 
mated to be of order ksTx = 0.003wQ. For compari- 
son we also plot the current in the non-interacting U = 
limit, calculated both exactly within the independent- 
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FIG. 3: Transport current J(t) for a quantum dot in the 
Kondo regime driven out of equilibrium by a step bias ap- 
plied to the left lead. J(t) is enhanced compared to the non- 
interacting (7 = case. The parameters used in the cal- 
culation are L = 100, t q = u;/(4v / 2), e q = —0.25 w, and 
U = oo, with a corresponding Kondo energy scale of order 
\zbTk = 0.003 w. The step in the bias is turned on at time 
t — 4t w within a rise time of t s = 2t w . The bias reduces 
the chemical potential of the left lead by a small amount, 
Sfi = —5 x 10~ 6 to. As the time advances small discrepan- 
cies show up between the more accurate M = 512 truncation 
(continuous curve with filled squares) and the M = 256 trun- 
cation (dotted curve). The units of current and time are the 
same as in Fig.jl]. 

particle picture, and via TdDMRG. The current for 
U = oo shows the expected Kondo enhancement p], p|. 
Upon decreasing the hopping amplitude between the dot 
and the leads to t q = 0.05w we find by contrast that the 
current at U — oo is less than at U = 0. We attribute 
this drop to the existence of a Coulomb blockade, which 
now overwhelms the vanishingly small Kondo resonance 
[ksTx — O(10~ 19 w)]. A more detailed discussion of our 
results, also in other settings, will be given elsewhere [Bp. 

In conclusion, we have shown that the DMRG algo- 
rithm can be extended to treat time-dependent problems. 
The TdDMRG method is a systematic tool to study the 
transient behavior of quantum many-body systems. The 
method has been illustrated by applying it to three sys- 
tems: spinless and spinful quantum dots and a junction 
between two Luttinger liquids. Our results agree with 
previous predictions but do not rely on uncontrolled ap- 
proximations. Reliable results can be obtained, even for 
systems driven far out-of-equilibrium. 
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